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Background:  Skin  color  changes  and 
mottling  are  frequently  described  signs  of 
hemorrhagic  shock  (HEM).  Based  on  this, 
we  developed  a  noninvasive,  noncontact  hy¬ 
perspectral  imaging  system  (HSI),  which 
quantifies  and  depicts  the  surface  tissue  sat¬ 
uration  of  oxygen  (SHS102)  for  each  pixel  in 
a  region  of  interest  (ROI).  Our  purpose  was 
to  assess  HSI  in  a  porcine  HEM  model.  We 
hypothesized  that  HEM  would  cause  de¬ 
creases  in  SHSI02  of  the  skin. 

IVlBthOdS:  The  HyperMed  HSI  sys¬ 
tem  employs  a  spectral  separator  to  vary 
the  wavelength  of  light  admitted  to  a  dig¬ 
ital  camera.  During  image  acquisition,  a 
“hypercube”  of  images,  each  at  a  separate 
wavelength,  is  generated  (at  5-nm  inter¬ 
vals,  from  500  to  600  nm).  Then,  the  visible 
light  spectrum  for  each  pixel  in  the  hypercube 
is  compared  by  linear  regression  to  standard 


spectra  for  oxyhemoglobin  (OxyHb)  and  de¬ 
oxyhemoglobin  (DeoxyHb).  The  resulting  fit 
coefficients  for  OxyHb  and  DeoxyHb  are 
used  to  calculate  SHSI02  values  for  each 
pixel  in  the  ROI.  The  mean  values  for 
OxyHb,  DeoxyHb,  and  SHSI02  across  the 
ROI  are  calculated.  Grayscale  SHS102  pic¬ 
tures  of  the  ROI  are  also  generated,  in 
which  the  brightness  of  each  pixel  is  pro¬ 
portional  to  its  value.  Seventeen  pigs, 
36.4  ±  0.11  kg,  underwent  standard  prep¬ 
aration,  and  were  maintained  on  ketamine 
and  isoflurane.  Normothermia  was  main¬ 
tained  (37°C  to  39°C).  The  hemorrhage 
group  (HEM,  n  =  9)  underwent  three 
blood  withdrawals,  each  10  mL/kg,  with 
15  minutes  between  withdrawals.  After 
the  third  withdrawal,  animals  were  resus¬ 
citated  with  lactated  Ringer’s  and  then 
shed  blood.  The  control  group  (CTRL, 


n  =  8)  received  intravenous  fluids  at  100 
mL/h.  HSI  images  were  obtained  of  the 
inner  hindlimb  throughout. 

Results:  All  HEM  animals  showed  lin¬ 
ear  decreases  in  both  mean  SHS,02  and 
OxyHb  values  with  blood  loss,  which  were 
reversed  by  resuscitation.  These  changes 
were  evident  on  the  grayscale  SHSI02  pic¬ 
tures,  but  not  to  the  naked  eye,  and  par¬ 
alleled  those  of  invasively  obtained  arte¬ 
rial  base  excess  and  mixed  venous  oxygen 
saturation. 

Conclusions:  HSI  is  a  promising 
noninvasive  and  noncontact  tool  for  quan¬ 
tifying  changes  in  skin  oxygenation  during 
HEM  and  resuscitation. 
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Changes  in  skin  color  and  mottling  of  the  skin  are  fre¬ 
quently  mentioned  as  important  clinical  findings  in 
critically  ill  patients:  “Decreased  perfusion  of  the  skin 
can  be  an  early  sign  of  shock .  .  .  Mottling,  pallor,  delayed 
capillary  refill,  and  peripheral  cyanosis  often  indicate  poor 
skin  perfusion.”1  However,  no  objective  method  of  quantify¬ 
ing  such  findings  currently  exists;  furthermore,  these  changes 
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may  be  absent  or  undetectable  to  the  naked  eye  in  a  large 
number  of  patients.  The  goal  of  this  study  was  to  assess  the 
utility  of  a  novel  technology,  hyperspectral  imaging  (HSI), 
for  noncontact,  noninvasive  monitoring  of  the  skin  during 
hemorrhagic  shock  and  resuscitation.  By  quantifying  and 
depicting  spatial  changes  in  hemoglobin  oxygen  saturation  in 
the  skin,  HSI  takes  advantage  of  the  fact  that  the  skin  is 
readily  available,  potentially  providing  a  “window”  through 
which  to  view  systemic  physiology.  This  application  of  HSI 
also  exploits  the  responsiveness  of  the  cutaneous  circulation 
to  hypovolemia. 

HSI  is  a  method  of  “imaging  spectroscopy.”  Spectros¬ 
copy  is  widely  used  to  monitor  metabolic  status  in  a  variety 
of  tissues.  For  example,  spectroscopic  methods  are  incorpo¬ 
rated  in  pulse  oximeters,  which  utilize  the  different  oxy-  and 
deoxyhemoglobin  absorption  bands  to  estimate  arterial  oxy¬ 
gen  saturation.2  Point  spectroscopy  in  the  near  infrared  range 
interrogates  the  tissue  saturation  of  oxygen  in  subcutaneous 
tissue,  muscle,  and/or  brain,  and  has  been  used  for  patient 
monitoring  in  hemorrhagic  shock,  in  extremity  compartment 
syndrome,  and  following  head  trauma.3-7 

Point  spectroscopy  systems  do  not  provide  information 
about  the  spatial  distribution  of  the  data.  Such  spatial  infor¬ 
mation  is  provided  by  HSI,  and  could  constitute  a  potential 
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Fig.  1.  Hyperspectral  imaging  generates  a  four-dimensional  data¬ 
set,  termed  an  image  cube  or  hypercube.  This  consists  of  two  spatial 
dimensions  (the  x  and  y  coordinates  of  each  pixel  in  the  region  of 
interest );  one  wavelength  dimension;  and  one  amplitude  dimension. 
In  other  words,  for  each  pixel  in  the  region,  the  visible  light  spec¬ 
trum  is  recorded. 

advantage  of  the  latter  approach.  In  HSI,  the  spectrum  of 
reflected  light  is  acquired  for  each  pixel  in  a  region,  and  each 
such  spectrum  is  subjected  to  standard  spectral  analysis  (Fig. 
1).  This  allows  the  creation  of  an  image  based  on  the  chem¬ 
ical  content  of  the  region  of  interest  (ROI). 

Used  for  decades  in  airborne  systems  for  the  analysis 
of  geological  features,  HSI  has  recently  been  applied  to 
biomedicine.8  It  has  been  employed  in  microscopic  studies  of 
histologic  sections.9'10  In  vivo,  HSI  has  been  used  to  demon¬ 
strate  the  macroscopic  distribution  of  skin  oxygen  saturation 
in  models  of  ischemia-reperfusion,  following  nitric  oxide 
inhalation  and  inhibition,  and  in  patients  with  sickle  cell 
disease  and  with  diabetic  foot  ulcers.11-14  One  application  of 
HSI  has  been  in  the  early  prediction  of  tissue  viability 
following  plastic  surgery:  tissue  that  has  insufficient  oxygen¬ 
ation  to  remain  viable  is  readily  apparent  from  oxygen  satu¬ 
ration  maps  calculated  from  near-infrared  spectral  images 
acquired  immediately  following  surgery.  By  contrast,  clinical 
signs  of  impending  necrosis  do  not  become  apparent  to  the 
naked  eye  for  6  to  12  hour  after  surgery.13,16  Assessment  of 
tissue  viability  following  burns  has  also  been  performed.17'18 

We  recently  developed  a  system  capable  of  hyperspectral 
imaging  in  the  visible  light  wavelength  range.  In  a  porcine 
pilot  study,  we  previously  used  this  system  to  observe 
changes  in  cutaneous  oxygen  saturation  following  chest 
trauma  and  hemorrhage,  which  were  not  evident  to  the  naked 
eye,  but  which  produced  hyperspectral  images  with  a  pro¬ 
nounced  mottling  pattern.19  Based  on  this  experience,  we 
hypothesized  that  the  image  intensity  of  hyperspectral  oxygen 
saturation  images  of  the  skin  would  decrease  during  hemor¬ 
rhagic  shock  in  a  porcine  model,  indicating  a  decrease  in 
oxygen  saturation  in  the  skin.  We  also  sought  to  gain  ini¬ 
tial  experience  with  HSI  during  hypothermia,  to  determine 


whether  this  common  complication  of  hemorrhagic  shock  and 
resuscitation  would  alter  the  HSI  images.  Our  ultimate  goal  is 
to  develop  a  system  that  demonstrates  blood-loss  induced 
changes  earlier  than  other  modalities,  or  which  provides  data 
which  are  otherwise  unavailable,  to  facilitate  rapid  diagnosis 
and  directed  treatment  of  shock. 

MATERIALS  AND  METHODS 

This  study  was  approved  by  the  Institutional  Animal 
Care  and  Use  Committee.  All  animals  were  cared  for  in 
accordance  with  the  guidelines  set  forth  by  the  Animal  Wel¬ 
fare  Act  and  other  federal  statutes  and  regulations  relating  to 
studies  involving  animals,  and  by  the  1996  Guide  for  the 
Care  and  Use  of  Laboratory  Animals  of  the  National  Re¬ 
search  Council.  All  animals  were  maintained  in  a  facility 
approved  by  the  International  Association  for  the  Assessment 
and  Accreditation  of  Laboratory  Animal  Care. 

Hemorrhagic  Shock  Study 

Seventeen  female  Yorkshire  pigs  (nine  in  the  hemorrhage 
group  and  eight  in  the  control  group),  weighing  36.4  +  0. 1 1  kg 
were  used  in  this  study.  The  animals  were  quarantined  for  1 
week  and  were  fasted  overnight  before  the  experiment. 

Animals  were  premedicated  with  250  mg  IM  Telazol 
(Lederle  Parenterals,  Inc.  Carolina,  PR).  After  induction  of 
anesthesia  with  isoflurane  delivered  through  a  mask,  they 
were  intubated,  and  were  placed  on  a  Datex-Ohmeda  anes¬ 
thesia  ventilator  with  a  tidal  volume  of  10  mL/kg  and  a 
respiratory  rate  of  12/min.  The  rate  was  adjusted  to  achieve 
normocapnea  (PaC02  =  35  to  45  mm  Hg).  Anesthesia  was 
maintained  with  a  mixture  of  isoflurane  (2%  to  2.5%)  and 
room  air.  Percutaneous  sheath  introducers  (8.5  Fr;  Arrow 
International,  Reading,  PA)  were  inserted  into  the  carotid 
artery  and  external  jugular  vein  bilaterally,  and  a  10-F  Foley 
catheter  was  inserted  into  the  urinary  bladder.  A  splenectomy 
was  performed  via  a  midline  laparotomy.  The  splenic  artery 
was  tied  off  before  splenectomy  to  allow  drainage  of  blood 
from  the  spleen  into  the  circulation.  An  infusion  of  lactated 
Ringer’s  solution  (LR)  at  1.5  times  the  spleen  weight  was 
administered  immediately  after  the  splenectomy.  At  the  end 
of  surgery  the  isoflurane  was  decreased  to  0.6%  and  an 
infusion  of  ketamine  (250  to  350  /rg/kg/min)  was  begun.  The 
ketamine-isoflurane  anesthesia  and  room-air  mechanical  ven¬ 
tilation  were  continued  until  the  end  of  the  study.  Depth- 
of-anesthesia  assessment  and  anesthetic  dose  adjustments 
were  made  as  needed.  Core  temperature  was  maintained 
between  37°C  to  39°C  by  means  of  an  external  heating  pad 
(Gayman  Industries,  New  York,  NY). 

A  flow-directed  pulmonary  artery  catheter  (VIP  Oximetry 
782HF75,  Edwards  Life  Sciences,  Irvine,  CA)  was  inserted  via 
the  external  jugular  vein  introducer  sheath  to  permit  measure¬ 
ment  of  mixed  venous  blood  gases  and  core  temperature.  One  of 
the  carotid  arterial  introducer  sheaths  was  used  for  measure¬ 
ment  of  the  arterial  blood  pressure  (ABP).  Clinical  pressure 
transducers  (Transpac  IV  trifurcated  monitoring  kit,  Abbott 
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Table  1  Experimental  Schedule,  Hemorrhage  Study 

Time  Point 

Event 

i 

Blood  transfusion 

2 

Bleed  1 

3 

Postbleed  1 

4 

Bleed  2 

5 

Postbleed  2 

6 

Bleed  3 

7 

Postbleed  3 

8 

LR  resuscitation 

9 

Post-LR 

10 

Blood  transfusion 

11 

Posttransfusion 

12 

End  of  study 

Bleed,  periods  during  which  1 0  ml/kg  blood  was  removed  (1 0  min 
each);  postbleed,  periods  of  observation  following  bleed  periods;  LR, 
lactated  ringer’s  solution;  blood  transfusion,  reinfusion  of  shed  blood. 


Critical  Care  Systems,  Abbott  Laboratories,  North  Chicago, 
IL)  were  used.  The  heart  rate  was  obtained  from  the  electro¬ 
cardiogram.  Regional  skin  temperature  was  monitored  on 
both  hindlimbs  using  thermocouples  (Physiotemp  BAT-  10, 
Type  T  thermocouple  thermometer  system,  Physiotemp  In¬ 
struments,  Clifton,  NJ). 

The  experimental  schedule  is  given  in  Table  1.  After  a 
postoperative  stabilization  period  (1  to  2  hours),  baseline  data 
were  collected.  The  animals  in  the  hemorrhage  group  then 
underwent  withdrawal  of  blood  through  the  carotid  line  with 
a  syringe.  Three  withdrawals,  each  10  mL/kg,  were  per¬ 
formed  at  a  constant  rate  of  1  mL/kg/min.  Blood  was  col¬ 
lected  into  a  bag  containing  CPDA  anticoagulant  (J-520 
blood  collection  unit,  Jorgensen  Laboratories,  Loveland, 
CO).  Each  of  the  three  10-minute  hemorrhage  periods  was 
followed  by  a  15-minute  observation  period.  Following  the 
third  observation  period,  the  animals  were  resuscitated  with 
intravenous  infusion  of  LR  at  1.5  times  the  shed  blood  vol¬ 
ume.  The  duration  of  the  LR  resuscitation  period  was  25  min¬ 
utes.  Additional  fluid  to  exceed  the  initial  resuscitation  volume 
was  then  administered  as  needed  to  return  the  heart  rate  and 
blood  pressure  toward  baseline  values.  The  25-minute  LR  infu¬ 
sion  period  was  followed  by  a  30-minute  observation  period, 
and  then  a  25-minute  period  during  which  the  shed  blood  was 
reinfused.  The  animals  were  observed  for  an  additional  hour 
thereafter  and  then  were  euthanized.  Blood  and  LR  were  in¬ 
fused  using  a  fluid  warmer  (Ranger  blood/fluid  warming 
system,  Model  254,  Augustine  Medical,  Prairie,  MN). 

Animals  in  the  control  group  underwent  similar  surgical 
preparation  and  received  a  maintenance  LR  infusion  at  100 
mL/h.  Data  were  obtained  at  the  following  time  points;  base¬ 
line,  after  each  10  mL/kg  blood  withdrawal,  after  LR  resus¬ 
citation,  and  after  blood  reinfusion. 

Hyperspectral  Imaging 

Hyperspectral  images  were  collected  using  a  10-bit 
Kodak  MegaPlus  1.6/  digital  camera,  consisting  of  a  1024  X 
1532  pixel  charge-coupled  device  (CCD)  element  and  a  10- 


bit  A/D  converter  (Roper  Scientific,  Trenton,  NJ).  The  cam¬ 
era  was  fitted  with  an  AF  50-mm  /72.8  EX  Macro  camera  lens 
(Sigma  Corporation,  Tokyo,  Japan).  Wavelength  selection 
was  provided  by  a  visible-wavelength,  liquid-crystal  tunable 
filter  (VIS2,  Cambridge  Research  and  Instrumentation, 
Woburn,  MA)  with  a  nominal  7-nm  bandpass.  This  filter  was 
fitted  to  the  front  of  the  camera  lens  (Fig.  2).  By  varying  the 
voltage  across  the  liquid  crystal  tunable  filter,  the  wavelength 
of  light  admitted  through  the  filter,  and  into  the  camera,  can 
be  varied.  A  100-W  quartz-halogen  light  powered  from  a 
Sorenson  150-W  DC  power  supply  (Elgar  Electronics,  San 
Diego,  CA)  was  used  to  provide  coaxial  illumination  for  the 
camera  system.  Custom  software  was  developed  under  the 
LabVIEW  system  (National  Instruments  Corporation,  Austin, 
TX)  to  control  the  image  acquisition  hardware  and  process. 

Images  of  the  proximal  inner  hindlimb  were  collected  at 
5-nm  intervals  from  500  to  600  nm,  with  the  integration  times 
adjusted  such  that  the  brightest  image  in  each  sequence  filled 
approximately  80%  of  the  full  well  capacity  of  the  CCD.  The 
data  resulting  from  one  of  these  sequences  is  known  as  an 
image  cube  or  “hypercube,”  each  of  which  requires  approx¬ 
imately  2  minutes  to  acquire  (Fig.  1).  Fiduciary  marks  were 
placed  on  the  skin  with  an  indelible  black  marker  to  facilitate 
image  registration. 

Images  were  collected  as  full  1024  X  1532  pixel  arrays,  and 
then  were  binned  down  to  256  X  382  arrays.  Background 
images  were  subtracted  from  each  image  in  the  image  cube.  The 
data  were  then  converted  to  optical  density  units  by  dividing  the 
sample  data  by  data  acquired  from  a  white  reflectance  standard 
using  a  Beer’s  Law  algorithm.  A  complete  description  of  this 
methodology  can  be  found  elsewhere20  and  in  Appendix  1. 

Imaging  System 


100  w 

quartz- 
halogen 
light 

Liquid 
Crystal 
Tunable 
Filter 
(LCTF) 

Fig.  2.  The  hyperspectral  imaging  system  consists  of  an  illumina¬ 
tion  source  and  a  digital  charge-coupled  device  (CCD)  camera.  In 
front  of  the  camera,  a  liquid  crystal  tunable  filter  is  used  to  vary  the 
wavelength  of  light  admitted  to  the  camera.  During  the  acquisition 
of  one  hypercube,  the  voltage  is  varied,  stepwise,  across  the  filter — a 
process  which  takes  about  2  minutes.  At  each  new  voltage,  an  image 
at  a  different  wavelength  of  light  is  obtained. 
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Briefly,  reference  oxyhemoglobin  and  deoxyhemoglobin  spectra 
were  obtained  in  electronic  format.21  A  four-term  linear  regres¬ 
sion  fitting  of  oxyhemoglobin,  deoxyhemoglobin,  offset  and 
slope  terms  was  then  performed  on  each  of  the  spectra  in  the 
image  cube.  The  regression  fit  coefficients  were  then  used  to 
calculate  the  relative  oxygen  saturation  percentage  for  each 
spectrum  in  the  image  cube: 

SHs[02  =  OxyHb/(OxyHb  +  DeoxyHb)  X  100, 

where  OxyHb  is  the  fit  coefficient  for  oxyhemoglobin,  and 
DeoxyHb  is  the  fit  coefficient  for  deoxyhemoglobin,  result¬ 
ing  from  the  linear  regression.19  The  mean  values  for  SHSI02, 
for  OxyHb,  and  for  DeoxyHb  were  calculated  across  the 
entire  ROI  of  each  image  at  each  timepoint,  yielding  vari¬ 
ables  which  were  then  subjected  to  further  statistical  analysis 
by  conventional  methods  (see  below).  In  addition,  grayscale 
images  were  generated  in  which  the  brightness  of  each  pixel 
is  proportionate  to  the  SHSI02  value  for  that  pixel.  All  spectral 
computations  were  done  using  Matlab  (The  Mathworks, 
Natick,  MA). 

Laser  Doppler  Imaging 

A  raster-optics-based  Laser  Doppler  Perfusion  Imager 
(Periscan  PIM  II,  Perimed  AB,  Stockholm,  Sweden)  was 
used  to  obtain  images  based  on  noncontact  laser  Doppler 
flowmetry  of  the  region  of  interest  (ROI).  Mean  perfusion 
was  calculated  off  line  for  16X16  pixel  windows  within  the 
ROI. 

Assays 

Arterial  and  mixed  venous  blood  samples  were  collected 
into  preheparinized  1-mL  disposable  syringes  at  each  sam¬ 
pling  time  point.  They  were  analyzed  using  an  AVL  Omni 
Modular  System  (Roche  Diagnostics  GmbH  Graz,  Austria). 

Statistical  Analysis 

Statistical  analysis  employed  SPSS  v.  10.1  (SPSS  Inc., 
Chicago,  IL)  and  Microsoft  Excel  2000  (Microsoft  Corp., 
Redmond,  WA).  Independent-sample  t  tests  were  used  to 
compare  control  versus  hemorrhage  animals  at  each  time- 
point.  Linear  regression  was  used  to  compare  continuous 
variables.  For  the  hyperspectral  variables,  bivariate  Pearson 
correlation  coefficients  were  calculated  for  each  animal  over 
time  in  the  hemorrhage  group,  from  baseline  to  time  point  7 
(postbleed  3).  Significance  was  accepted  at  p  <  0.05.  Data 
are  presented  as  means  ±  SEM. 

RESULTS 

As  can  be  seen  in  Table  2,  blood  withdrawal  resulted  in 
an  early  drop  in  systolic  arterial  pressure,  which  became 
statistically  significant  in  comparison  with  the  control  group 
during  the  first  withdrawal,  and  which  remained  decreased 
until  after  LR  resuscitation.  The  heart  rate  increased  with 
blood  withdrawal,  although  in  delayed  fashion — becoming 


significant  only  during  the  second  bleed;  it  returned  to  control 
levels  after  reinfusion  of  shed  blood.  Other  physiologic  pa¬ 
rameters  are  given  in  Table  2. 

Oxygen  saturation  images  in  which  the  brightness  of 
each  pixel  is  proportionate  to  the  intensity  of  the  SHSI02  for 
that  pixel  are  shown  in  Figure  3.  Here,  both  baseline  images 
and  images  obtained  during  the  third  postbleed  period  are 
included.  The  decrease  in  image  intensity  is  evident  on  these 
images.  Qualitatively,  some  animals,  but  not  all,  demon¬ 
strated  an  increase  in  mottling  during  shock,  also  evident  on 
the  oxygen  saturation  images.  Neither  of  these  changes  was 
evident  to  the  naked  eye.  Quantitatively,  the  mean  SHSI02 
value  for  these  images  decreased  linearly  with  blood  with¬ 
drawal,  becoming  significantly  decreased  in  comparison  with 
control  animals  after  the  third  bleed,  and  restored  to  control 
levels  by  resuscitation  (Fig.  4,  Table  2).  The  mean  value  of 
the  OxyHb  fit  coefficient  for  the  ROI  also  decreased  linearly 
with  hemorrhage,  but  it  showed  an  earlier  statistically  signif¬ 
icant  decrease,  after  the  second  bleed,  which  was  also  re¬ 
stored  by  resuscitation  (Fig.  5,  Table  2).  Meanwhile,  the 
mean  DeoxyHb  fit  coefficient  for  the  ROI  trended  upward 
during  and  after  the  third  bleed,  but  these  changes  were  not 
significant  (Fig.  6,  Table  2). 

Pearson  correlation  coefficients  were  calculated  for  each 
animal’s  individual  hyperspectral  variables  over  time  points  1 
to  7  (baseline  to  postbleed  3).  These  results  are  given  in  Table 
3.  Note  that  the  correlation  coefficients  for  the  primary  end¬ 
point  variable,  SHSI02,  ranged  from  —0.755  to  —0.976  (mean 
—0.863);  in  other  words,  this  variable  showed  the  expected 
linear  decrease  in  SHSI02  with  hemorrhage  in  all  cases. 
Likewise,  all  but  one  animal  demonstrated  the  expected  linear 
decrease  in  OxyHb  with  hemorrhage.  By  contrast,  results  for 
DeoxyHb,  although  positive  as  expected  for  all  animals,  were 
more  variable. 

The  SD  of  the  SHSI02  images,  a  rough  index  of  the 
degree  of  mottling  present,  did  not  change  significantly  (data 
not  shown).  Laser  Doppler  imaging  demonstrated  a  linear 
decrease  in  skin  blood  flow  for  the  ROI  with  blood  loss, 
which  became  significant  after  the  third  bleed  and  which  was 
reversed  by  resuscitation  (Table  2). 

Linear  regression  was  performed  to  examine  the  possible 
relationship  between  mixed  venous  saturation  of  oxygen  and 
SHsi02-  This  analysis  demonstrated  a  linear  relationship,  with 
a  relatively  low  r2  of  0.12  (p  <  0.001,  df  =  1 14).  Correlation 
with  the  systolic  arterial  pressure  was  similar  (r2  =  0.14,  p  < 
0.001,  df  =  202).  Although  the  laser  Doppler  image  intensity 
appeared  to  follow  a  time  course  similar  to  that  of  SHSI02, 
there  was  no  relationship  between  the  two  variables  on  linear 
regression  (r2  =  0.01,/?  =  0.312,  df  =  112). 

DISCUSSION 

The  principal  finding  in  this  study  is  that  the  mean 
intensities  of  both  SHSI02  and  OxyHb  images  of  the  skin, 
obtained  by  hyperspectral  imaging  in  the  visible  wavelength 
range,  decreased  during  hemorrhagic  shock  and  were  restored 
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Hyperspectral  Imaging 
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Baseline  Shock  Baseline  Shock 


Fig.  3.  These  images  were  generated  by  calculating  the  SHSI02  value  for  each  pixel  in  the  field  of  view.  Then,  each  pixel  was  given  a 
grayscale  brightness  value  proportionate  to  the  SHSI02  value.  The  resulting  images  are  shown.  Clear  differences  in  SHSI02  are  evident 
between  baseline  images  and  those  produced  after  the  third  bleed.  In  addition,  some  animals  demonstrated  an  increase  in  regional 
heterogeneity,  or  mottling — most  notably  the  third  animal  on  the  left.  These  changes  were  not  evident  to  the  naked  eye. 


during  resuscitation  in  anesthetized  pigs.  These  changes 
roughly  paralleled  those  observed  in  several  invasively  obtained 
variables,  to  include  the  systolic  arterial  pressure,  the  arterial 
base  excess,  and  the  mixed  venous  saturation  of  oxygen.  Thus, 
HSI  may  be  a  useful,  noninvasive  method  of  monitoring  the 
response  to  shock  and  the  adequacy  of  resuscitation. 


SHSI02  (mean) 


Timepoint 

1  I  I  I  I 

Bleed  1  Bleed  2  Bleed  3  LR  Blood 

Fig.  4.  The  mean  SHSI02  for  the  region  of  interest  was  calculated 
and  is  shown  here.  SHSI02  decreased  linearly  with  blood  withdrawal 
and  became  significantly  different  from  control  values  after  the  third 
bleed. 


Systolic  blood  pressure  dropped  more  rapidly  than  did 
the  other  variables,  including  the  mean  SHSI02  and  OxyHb 
variables.  The  rapidity  of  this  drop  in  blood  pressure  during 
relatively  moderate  blood  loss,  coupled  by  a  delay  in  the 
compensatory  tachycardic  response  to  the  decrease,  indicates 
that  cardiovascular  compensation  for  hemorrhage  may  have 
been  impaired  in  this  anesthetized  porcine  model.  The  de¬ 
pressive  effects  of  general  anesthesia  on  sympathetic  nervous 

HSI  Oxyhemoglobin  Coefficient 


Timepoint 

I  i  i  i  i 

Bleed  1  Bleed  2  Bleed  3  LR  Blood 

Fig.  5.  The  mean  OxyHb  fit  coefficient  for  the  same  region  of 
interest  is  shown  here.  This  value  also  decreased  linearly  with 
hemorrhage  and  became  significantly  different  from  control  values 
after  the  second  bleed. 
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HSI  Deoxyhemoglobin  Coefficient 


i  i  i  i  i 


Bleed  1  Bleed  2  Bleed  3  LR  Blood 

Fig.  6.  The  mean  DeoxyHb  fit  coefficient  increased  with  hemor¬ 
rhage,  but  this  change  was  not  statistically  significant  in  comparison 
with  control  animals. 

system  activity,  on  peripheral  vascular  responsiveness,  on 
cutaneous  vasoconstriction,  on  baroreflex  sensitivity,  and  on 
cardiac  function  are  well  known. 22-26  These  effects  limit  the 
applicability  of  models  which  employ  general  anesthesia, 
such  as  that  used  in  the  present  study,  to  conscious  human 
trauma  patients.  Recently,  we  have  improved  the  image  reg¬ 
istration  capabilities  of  the  HSI  system,  such  that  immobili¬ 
zation  of  subjects  by  means  of  general  anesthesia  should  not 
be  required  for  future  studies.  Whether  studies  in  unanesthe¬ 
tized  subjects  would  reveal  increased  sensitivity  to  lower 
volumes  of  blood  loss,  and  whether  HSI  would  be  superior  to 
traditional  vital  signs  under  those  circumstances,  remains  to 
be  determined. 

Mottling  was  inconsistently  observed  in  these  hyperspec¬ 
tral  images,  in  contrast  to  the  pronounced  increase  in  mottling 
which  we  previously  observed  in  a  single  case  of  chest 
trauma  and  hemorrhage.19  The  reasons  for  this  are  unclear. 


Table  3  Pearson  Correlation  Coefficients  and 
Hyperspectral  Data 


Subject 

OxyHb 

PC 

Mean 

Sig 

DeoxyHb 

PC 

Mean 

Sig 

ShS|02 

PC 

Mean 

Sig 

i 

-0.898 

0.006 

0.659 

0.108 

-0.872 

0.011 

2 

0.562 

0.189 

0.874 

0.010 

-0.755 

0.050 

3 

-0.994 

<0.001 

0.107 

0.819 

-0.970 

<0.001 

4 

-0.727 

0.064 

0.792 

0.034 

-0.828 

0.021 

5 

-0.954 

0.001 

0.648 

0.115 

-0.921 

0.003 

6 

-0.653 

0.112 

0.557 

0.194 

-0.704 

0.077 

7 

-0.922 

0.003 

0.925 

0.003 

-0.939 

0.002 

8 

-0.976 

<0.001 

0.401 

0.373 

-0.976 

<0.001 

9 

-0.905 

0.005 

0.077 

0.870 

-0.804 

0.029 

For  each  animal  in  the  hemorrhage  group  were  analyzed  over 
time  points  1  to  7. 

OxyHb  mean,  mean  value  of  the  oxyhemoglobin  fit  coefficient  for 
the  region  of  interest;  DeoxyHb  mean,  mean  value  of  the  deoxyhe¬ 
moglobin  fit  coefficient  for  the  region  of  interest;  SHS|02  mean,  mean 
value  of  the  oxygen  saturation  of  hemoglobin  for  the  region  of  interest; 
PC,  Pearson  correlation  coefficient;  sig,  two-tailed  significance  level. 


However,  anesthesia  may  also  impair  the  control  mechanisms 
for  cutaneous  blood  flow,  causing  a  dampening  of  the  blood 
flow  oscillations  which  normally  occur  in  response  to  hypo¬ 
volemia.  Braverman  and  colleagues  described  the  relation¬ 
ship  between  the  microanatomy  of  the  cutaneous  circulation 
and  regional  heterogeneity  in  blood  flow  by  laser  Doppler 
flowmetry. 27-30  They  also  described  synchronicity  in  cutane¬ 
ous  blood-flow  oscillations  between  sites  on  ipsilateral  and 
contralateral  limbs,  suggesting  that  such  oscillations  are  con¬ 
trolled  centrally  by  the  sympathetic  nervous  system.  Indeed, 
stimulation  of  the  lumbar  sympathetic  trunk  at  appropriate 
frequencies  increases  the  amplitude  of  skin  blood-flow  oscil¬ 
lations  measured  by  laser  Doppler.  ’1  Such  increases  in  vaso- 
motion  may  play  a  protective  role,  by  enabling  red  blood  cells 
to  traverse  the  microvasculature  under  conditions  of  in¬ 
creased  vascular  resistance.32  Further  studies  employing  HSI 
in  unanesthetized  humans  undergoing  lower  body  negative 
pressure  are  being  conducted,  and  may  determine  whether  an 
increase  in  regional  heterogeneity  in  SHSI02  images  is  a 
prominent  feature  of  the  response  to  central  hypovolemia.  In 
addition,  novel  image  analysis  methods  of  detecting  and 
quantifying  such  mottling  are  being  developed.  Finally,  direct 
measurement  of  cutaneous  sympathetic  nerve  activity  would 
shed  light  on  control  of  the  cutaneous  circulation  during 
shock.33 

Thermoregulation  is  a  major  function  of  the  cutaneous 
circulation,34’35  and  hypothermia  is  commonly  present  in 
trauma  casualties.  Thus,  the  relationship  between  systemic 
hypothermia  and  localized  cutaneous  oxygen  saturation  is 
relevant  to  an  assessment  of  HSI  technology.  We  were  con¬ 
cerned  that  hypothermia  would  activate  cutaneous  sympa¬ 
thetic  nerves  and  cause  a  reduction  in  cutaneous  blood  flow; 
in  fact,  this  response  can  be  used  as  a  test  of  the  integrity  of 
the  autonomic  nervous  system.36’37  In  turn,  a  reduction  in 
cutaneous  blood  flow  could  cause  a  decrease  in  skin  oxygen 
saturation,  assuming  local  metabolic  demands  remained  con¬ 
stant.  The  decrease  in  mean  SHSI02  and  OxyHb  values  during 
the  hemorrhagic  shock  study  were  not  caused  by  systemic 
hypothermia;  core  and  skin  temperature  were  maintained 
during  the  experiment. 

In  a  separate  pilot  study,  moderate  and  profound  hypo¬ 
thermia  did  not  appear  to  produce  similar  effects  (data  not 
shown).  Possibly,  decreased  metabolic  demands  defended 
against  desaturation  during  moderate  and  profound  hypother¬ 
mia.  However,  further  research  would  be  needed  to  evaluate 
fully  the  effect  of  hypothermia  on  HSI. 

The  technology  used  in  this  study  performed  imaging  in 
the  visible  wavelength  range  of  500  to  600  nm.  This  wave¬ 
length  selection  provides  improved  signal-to-noise  ratio  over 
near-infrared  spectroscopy  for  the  measurement  of  hemoglo¬ 
bin  spectra,  because  the  hemoglobin  absorption  signal  is 
much  stronger  in  the  visible  than  in  the  near-infrared  range. 
The  utility  of  near-infrared  point  spectroscopy  lies  in  its 
ability  to  provide  information  about  subsurface  tissue  oxy¬ 
genation,  for  example  in  brain  or  muscle,  and  which  derives 
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from  the  increased  tissue  penetration  of  near-infrared  light !  7 
By  contrast,  visible  light  spectroscopy,  as  employed  in  this 
study,  interrogates  hemoglobin  saturation  at  a  more  superfi¬ 
cial  level,  most  likely  within  dermal  capillaries.  Near-infrared 
spectroscopy  can  also  be  used  to  measure  tissue  water 
content.38  Thus,  visible  and  near  infrared  approaches  provide 
potentially  complementary  information.  Furthermore,  tech¬ 
nology  that  simultaneously  performs  hyperspectral  imaging 
in  both  the  near  infrared  and  visible  ranges  is  under  devel¬ 
opment.  In  darkly  pigmented  subjects,  epidermal  melanin 
dominates  the  visible  light  spectrum.  Several  studies  suggest 
that  this  may  not  impair  the  ability  of  HSI  systems  to  char¬ 
acterize  changes  in  cutaneous  oxygen  saturation.1214  39 

Calibration  of  this  HSI  system  was  performed  using  refer¬ 
ence  spectra.21  Direct  measurement  of  OxyHb  and  DeoxyHb 
levels  in  skin  has  not  been  performed  in  conjunction  with  HSI. 
It  can  be  seen  (Figs.  5  and  6)  that  there  was  more  variability  in 
DeoxyHb  than  in  OxyHb  levels  in  control  animals,  and  in  both 
values  in  control  animals  as  compared  with  hemorrhaged  ani¬ 
mals.  Whether  the  changes  in  OxyHb  and  DeoxyHb  levels  seen 
in  this  study  are  clinically  significant — whether,  for  example, 
certain  values  are  associated  with  survival — can  only  be  deter¬ 
mined  by  more  work  in  animals  and  humans. 

Laser  Doppler  imaging  is  another  approach  to  imaging 
the  skin  during  shock  and  resuscitation.  Light  and  colleagues 
have  envisioned  the  use  of  this  technology  as  a  guide  to 
resuscitation  following  thermal  injury.40  HSI  and  laser  Dopp¬ 
ler  imaging  are  different  technologies,  the  former  providing 
information  about  tissue  oxygenation  and  the  latter  about 
blood  flow.  Further  improvements  in  HSI  may  offer  advan¬ 
tages  over  currently  available  laser  Doppler  imagers  with 
respect  to  image  resolution  and  acquisition  time.  In  conclu¬ 
sion,  we  have  developed  a  hyperspectral  imaging  system 
which  demonstrates  changes  in  skin  oxygen  saturation  during 
hemorrhagic  shock  and  resuscitation.  Further  studies  are  planned 
to  elucidate  fully  the  potential  of  this  new  technology. 

APPENDIX  1:  IMAGE  PROCESSING 

The  following  is  a  more  detailed  description  of  the  HSI 
process  used  in  this  study,  adopted  from  Gillies  et  al.19  Raw 
sample  images  (S)  and  reference  images  (R)  were  smoothed 
using  a  three-point  mean  filter  along  spectral  dimension  of 
the  image  cubes.  The  smoothed  sample  image  cubes  (S)  were 
ratioed  against  the  smoothed  reference  image  cubes  (R)  to 
convert  the  data  to  an  optical  density  scale  (A).  The  optical 
density  at  pixel  i,j  (Aij)  of  each  raw  image  was  calculated  by: 

Ay  =  “  lOgl()  (Sy/Ry)  (l) 

The  conversion  of  the  data  to  optical  density  allowed  the 
spectra  to  be  analyzed  using  standard  spectroscopic  algo¬ 
rithms 20  Linear  regression  fit  coefficients  cl,  c2,  and  c3  were 
calculated  for  each  spectrum  in  an  image  cube  from  in  vivo 


reference  spectra21  for  oxyhemoglobin,  deoxyhemoglobin, 
and  an  offset  spectrum,  respectively,  by: 

Sy  =  \\ciOxyHb  +  c2DeoxyHb  +  c2Offset\\2  (2) 

These  coefficients  allowed  the  determination  of  total 
hemoglobin  concentration  (THb): 

Total  Hemoglobin  =  THb:j  =  coxyij  +  cdeom  (3) 
and  relative  percent  saturated  oxygen  (02sat): 

%02patv  =  [coxyJ ( THbjj)]  X  100  (4) 

at  each  pixel  location.  Grayscale  images  based  on  these  values 
were  analyzed  using  standard  image  analysis  techniques.  SHSI02 
is  calculated  as  the  mean  value  of  %02sat  across  the  region  of 
interest. 
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